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We argue that an arbitrary general relativistic static anisotropic fluid sphere, (static and 
spherically symmetric but with transverse pressure not equal to radial pressure), can 
nevertheless be successfully mimicked by suitable linear combinations of theoretically 
attractive and quite simple classical matter: a classical (charged) isotropic perfect fluid, 
a classical electromagnetic field, and a classical (minimally coupled) scalar field. While 
the most general decomposition is not unique, a preferred minimal decomposition can 
be constructed that is unique. We show how the classical energy conditions for the 
anisotropic fluid sphere can be related to energy conditions for the isotropic perfect 
fluid, electromagnetic field, and scalar field components of the model. Furthermore we 
show how this decomposition relates to the distribution of both electric charge density 
and scalar charge density throughout the model. The generalized TOV equation implies 
that the perfect fluid component in this model is automatically in internal equilibrium, 
with pressure forces, electric forces, and scalar forces balancing the gravitational pseudo¬ 
force. Consequently, we can build theoretically attractive matter models that can be used 
to mimic almost any static spherically symmetric spacetime. 
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1. Introduction 

Perfect fluid spheres in general relativity have almost a full century of history, dating 
back to Schwarzschild’s interior solution of 19161^ (corresponding to a spatially 
uniform density with a spatially varying but isotropic pressure). More recently, some 
particularly important contributions are the Delgaty-Lake review article of 1998P 
and the subsequent rapid (and extensive) development of “algorithmic techniques” 
in the first decade of the 21 st century PEH 

However, when it comes to dealing with anisotropic fluid spheres, (where the 
radial pressure need not equal the transverse pressure), the situation is considerably 
messier. See, for instance, references fl5H29l Such anisotropic fluid spheres are some¬ 
times believed to be physically relevant for the description of neutron stars, (see, 
for example, references i30lf34l) . and other compact objects, and are known to be 
relevant for the internal structure of gravastars E5H3H At a pinch the formalism can 
also readily be adapted to anisotropic solid spheres, such as idealized spherically 
symmetric planets, or the possibly solid crust of a neutron star. Typically one can 
proceed only at the high cost of making some very specific ansatz; either for some 
specific metric component, or for some specific relationship between the compo¬ 
nents of the stress-energy tensor. Often these ansatze are less than well physically 
motivated. 

In the current article we shall endeavour instead to build a reasonably general 
theoretical model to “mimic” a generic anisotropic fluid sphere. We shall do this 
by “mimicking” the stress-energy of a generic anisotropic fluid sphere in terms of 
some linear combination of “charged perfect fluid” plus “electromagnetic field” plus 
“massless scalar field”. By construction, the perfect fluid component in this model 
is automatically in internal equilibrium, with pressure (buoyancy) forces, electric 
forces, and scalar forces exactly balancing the gravitational pseudo-force. 

We shall investigate the extent to which such a decomposition can be carried 
out, the extent to which such a decomposition is (or can be made to be) unique, 
and the extent to which we can say something concerning the spacetime geometry. 
We shall also investigate the extent to which the classical energy conditions®®!! 
are (or can be made to be) satisfied, and the extent to which we can pin down the 
(effective) electric charge density and scalar charge density that are used to “mimic” 
the anisotropies in the stress-energy. Overall, we will seek to build a simple and 
straightforward theoretical model that can be usefully used to “mimic” a very large 
class of anisotropic fluid spheres. 

For definiteness we shall adopt the conventions that the spacetime metric is 
presented in the usual Schwarzschild curvature coordinates 

r\r 2 

ds 2 = -e~ 2 ^ r) dt 2 + - ——- + r 2 ( dd 2 + sin 2 9 d^ 2 ), 

1 — 2 m[r)/r 

and that in the associated orthonormal basis the total stress-energy is: 


(i) 
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2. Stress-energy tensors 

The relevant stress-energy tensors we use are utterly standard, see for instance 
references [3D and [45] 


• For a perfect fluid we have 

T { ab = (pi+P{)V a V b +p { g ab . 

• For the electromagnetic field we have 


jiab 
^ em 


= F ac g cd F bd - ±g ab (F cd F cd ). 


• For the (minimally coupled) massless scalar field we have 

~ab ( „cd. „ 


T s ab = - ^g ab (g cd Md)- 


(3) 

(4) 

(5) 


rj-iab _ 


( 6 ) 


F hb = E 


(7) 


Once one restricts to spherical symmetry, and adopts an orthonormal basis, the 
stress-energy for a perfect fluid is: 

Pi 0 0 0 
0 p { 0 0 
0 0 pi 0 

0 0 Oft 

Similarly, spherical symmetry plus an orthonormal basis imply that the electromag¬ 
netic field-strength tensor is: 

0+100 
-1 0 0 0 
0 0 0 0 
0 0 0 0 

Consequently, spherical symmetry, an orthonormal basis, and tracelessness implies 
that the electromagnetic stress-energy is: 

+10 0 0 
0-100 
0 0+10 
0 0 0 +1 

Finally, using spherical symmetry and an orthonormal basis, for a scalar field: 

0 0 0 0 

</> ; V = (Vc/>) 2 


rpab _ 

1 em 


l 


-E 1 


( 8 ) 


0+100 
0 0 0 0 
.0 0 0 0 

Consequently, the stress-energy for a massless minimally coupled scalar is: 

'+1000 


(9) 


T s ab = ^(Vcf >) 2 


0+100 
0 0-10 
0 0 0 -1 


( 10 ) 
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Note that for the electromagnetic field we automatically have (p r ) em < (pt)em, 
while in contrast for the scalar field (p r ) 8 > (pt)s- It is this relative magnitude flip 


that is essential in mimicking an arbitrary anisotropic fluid sphere. 

Combining these results we now have: 

p = p i+ l -E 2 + l -(S ( j>) 2 -, 

(11) 

Pr = Pi ^ E 2 + i(V0) 2 ; 

(12) 

Pt =Pf + ^E 2 - i(V0) 2 . 

(13) 

Then in particular 

Pr+Pt= 2pf; and p r —p t = (V^) 2 - E 2 . 

(14) 

One specific way of inverting this, (by far the simplest), is to take: 

Pi = \{Pr+Pt ); 

(15) 

Pi = P- \ \Pr Pt |; 

(16) 

(V</>) 2 = ma x{p r -p t , 0}; 

(17) 

E 2 = ma x{p t - p r , 0}. 

(18) 


By using this decomposition we see that any static and spherically symmetric stress 
tensor can be mimicked by a linear superposition of 

(perfect fluid) + (electromagnetic) + (massless scalar). 

Indeed: 

• If p r < pt everywhere within the anisotropic fluid sphere, then one can get 
away with only using (perfect fluid) + (electromagnetic) contributions. 

• If p r > Pt everywhere within the anisotropic fluid sphere, then one can get 
away with only using (perfect fluid) + (massless scalar) contributions. 

Note that in regions where the transverse pressure is greater than the radial pressure 
the model has |V</>| = 0, so the scalar field is constant. (We shall soon see that the 
scalar charge density is zero in such regions.) In such regions E(r) ^ 0, and the 
net charge inside the sphere of radial coordinate r is simply Q(r) = AnE(r)r 2 . 
Unless E(r) is fine-tuned so that E(r) oc 1/r 2 , there will be a non-zero electric 
charge density in such regions. The electrically charged perfect fluid is then subject 
to both electric and pressure (buoyancy) forces, as well as the gravitational pseudo¬ 
force. Because the model is (by construction) static, these forces will automatically 
balance each other so that the model is in internal equilibrium — more details on 
this point below, when we discuss the generalized TOV equations. 



October 2, 2015 0:46 anisotropic-spheres-2015-october-01 


Mimicking static anisotropic fluid spheres in general relativity 5 

If the surface of the anisotropic fluid occurs in such a region, where the transverse 
pressure is greater than the radial pressure, then in order to match an exterior 
Schwarzschild geometry the electric field E(r) would have to be discontinuous at 
the surface. (That is, E(r~) ^ 0, but E(rf ) = 0.) This implies the need for a surface 
charge density to cancel the immediate sub-surface electric field. This is merely a 
slightly unusual feature of our model, it is in no sense problematic. 

If one violently objects to the presence of surface charge density, then one could 
instead simply match the anisotropic sphere to a Reissner-Nordstrom spacetime; 
however that would leave the anisotropic fluid sphere with a net electric charge, 
which we would expect to be rapidly neutralised by quite standard astrophysical 
processes. Overall, it seems more appropriate to work with a non-zero surface charge 
density but net overall charge of zero. 

In counterpoint, in those regions where the radial pressure is greater than the 
transverse pressure our model has E = 0; this corresponds to both zero electric 
charge density and the somewhat stronger requirement that the net electric charge 
(interior to this region) be zero. Since V<^ ^ 0 in such regions there will now typically 
be a non-zero scalar charge density. The scalar-charged perfect fluid is then subject 
to both scalar forces and pressure (buoyancy) forces, as well as the gravitational 
pseudo-force. Because the model is (by construction) static, these forces will auto¬ 
matically balance each other — more details on this point below. If the surface of the 
anisotropic fluid occurs in such a region then to match an exterior Schwarzschild ge¬ 
ometry the scalar field </>(r) must be continuous at the surface, and constant outside 
the anisotropic fluid sphere. The matching to an exterior Schwarzschild geometry 
is in this case trivial. 

3. Stress-energy tensors and their covariant divergences 

We now proceed to set up a generalized TOV system of equations, ultimately based, 
(as usual), on the covariant conservation T ab . b = 0 of the total stress-energy. (This 
allows us to study the internal forces acting on the various matter components in 
our model, and so directly address questions of internal equilibrium.) To do so, let 
us first calculate the covariant divergence for each component separately. 

For a perfect fluid it is a standard result that®®] 

T { ab , b = {( Pi + P i)V a ;b V b + (g ab +V a V b )(p { ). b } + {( P V b ), b + P (V b , b )} V a . (19) 

Here the two terms are by construction 4-orthogonal. One of these terms relates 
the 4-acceleration (in the present context, the gravitational pseudo-force) to the 
pressure gradient (buoyancy forces), the other ultimately provides a conservation 
law. 

For the electro-magnetic field one has 

T e t, b = F ac g cd {F bd , b ) + g am F md , b F bd - \g ab F cd , b F cd . 


( 20 ) 
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Use of the Bianchi identity for F, and a few standard manipulations, quickly reduces 
this toim 

T e t ;b = - F ac g cd (F db ;b )• (21) 

Now writing this in terms of the 4-current 

(F db , b ) = J d = a em V d , (22) 

where a em is the electric charge density, we have the standard Lorentz force law^HUl 

T e t, b = -F ab (a em v b ) = —(Jem F ab V b . (23) 

Note that, simply because we have assumed everything is static, the 4-current is 
automatically parallel to the 4-velocity of the perfect fluid component. 

For the (minimally coupled) massless scalar field a brief computation yields®!® 

T s ab ;b = <p a .,b<P b + <t> ia <j>’ b -,b - g ab (g cd Md) = 0 ;a (0 ;b ;ft)- (24) 

Now set {4>' b -, b ) = (J s , where <j s is the scalar charge density. Then the scalar force 
density is given by 

T s ab . b = a s ft a . (25) 

Finally, combine the three (perfect fluid), (electro-magnetic field), and (scalar field) 
components. In the absence of other forms of matter, because the total stress-energy 
must be covariantly conserved, T ab . b = 0, we have 

{(p f +Pt)V a . b V b + ( g ab + V a V b )( Pi ), b } + {(p f V b ), b + Pl (V b , b )} 

~F ab (a em V b ) + a s <t>’ a = 0. (26) 

Projecting this first along the fluid 4-velocity V, and then perpendicular to V, we 
have the two equations 

(p f V b ) ;b + P! (V b . b ) = -a s ct>, a V a , (27) 

and 

(Pf + Pf)V a ;b V b + ( g ab + V a V b )([p f ] ;b + a s cf . b ) - F ab (a em V b ) = 0. (28) 

These two equations nicely summarize how the scalar charge and electric charge 
interact with the perfect fluid pressure, perfect fluid density, and 4-acceleration. 
In particular, recall that the spacetime (and the matter content) has been assumed 
static, (so the timelike Killing vector K a is parallel to the fluid 4-velocity V a ). Then 
things simplify nicely. First of all, < j)- a V a = 0, but we also have both p- a V a = 0 and 
p, a V a = 0. This then leads to V a - a = 0, and in fact the first equation (l27l) is vacuous. 
However the second equation (1281) is very definitely non-vacuous — it simplifies only 
slightly to yield 

(Pf + Pi)V a , b V b + g ab ([p f] ;6 + a s cf, b ) - F ab (a em V b ) = 0. 

This will quickly lead to the generalized TOY equations. 


(29) 
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4. Generalized TOV system of equations 

Now consider again equation (l29l) . and note that due to spherical symmetry the 
three individual pieces of this equation can point only in the radial direction. Let 

n a = (0,y/l - 2m(r)/r,0,o) (30) 

be the unit vector in the radial direction. Then from equation (1291) we have 

- (pf + Pi)V a V b n a -b + n“([Pf];6 + cr s (j)-b) + (JemE = 0. (31) 


But we know that if tr s = 0 = a em this must be equivalent to the usual TOV system 
of equations 


dpf _ (pf+pf)(m + 47rp f r 3 ) 
dr r 2 (l — 2 m/r) 


dm 

dr 


47rp p r 2 . 


(32) 


Consequently, reinstating the charge densities by taking er s ^ 0 ^ cr em , we see 


dpf n d</> _ (Pf +Pf)(w + 47rp f r 3 ) _ a em E 
dr s dr r 2 (l-2m/r) ^/l - 2m/r' 


(33) 


That is, the generalized TOV equation is now 

dpf (pf +Pf)(m + 47rp f r 3 ) cr em E d </> 

dr r 2 (l — 2m/r) y/Y^2mJr ^ dr ’ 


(34) 


supplemented by the generalized mass relation 


d?n 

dr 


= 4npr 2 = 4n(p f + p em + Ps)r 2 . 


(35) 


Note that the first of these two equations explicitly and automatically balances 
the pressure force, electric force, and scalar force against the pseudo-force due to 
gravity. (This is as it must be, since we started from a static model there can be no 
motion of the perfect fluid component.) 


Since this is a tricky point let us be explicit: Simply by construction the model 
is automatically in internal equilibrium. We cannot however, as yet say whether 
this equilibrium is stable, neutral, or unstable — at least not without additional 
hypotheses, (such as some equation of state relating pressure, the mass density, and 
the charge densities). In general, the perfect fluid need not, of course, be homoge¬ 
neous, but this was already true in the Schwarzsclrild interior geometry of 1916, and 
so should be no particular surprise in this more general context. Similarly, while the 
Schwarzschild interior solution is by construction automatically in internal equilib¬ 
rium, we cannot say whether this equilibrium is stable, neutral, or unstable without 
some additional hypothesis. In this regard the model we are adumbrating in the 
current article is certainly no worse than the usual models for perfect fluid spheres. 
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This generalized system of TOV equations still has to be supplemented with the 
equations of motion for the electromagnetic and scalar fields. Specifically, for the 
electric field we have (in the quite usual fashion) 


Here Q(r) is the total net charge interior to radius r. Explicitly 



(37) 


For the scalar field we have A <f) = a Sl which reduces to 



(38) 


That is 



(39) 


This can now be integrated to explicitly yield 



(40) 


If one desires an explicit formula for the scalar field itself then, up to an arbitrary 
additive constant, we have 



(41) 


This is now a complete set of differential and integral equations which uniquely 
determine the three-component model used to mimic the anisotropic fluid sphere. 

5. Classical energy conditions 

Let us now consider the classical energy conditions, see for instance references I37W431 
Since the electromagnetic and scalar field stress energy tensors both satisfy all the 
classical energy conditions, it is clear that provided the fluid component satisfies 
the classical energy conditions, then the total stress-energy must also satisfy the 
energy conditions. (The converse need not always hold.) 

Specifically, let us consider the NEC and note that 



(42) 


This implies either 


Pr> Pt- 


Pf+Pf = p+pt', 


(43) 


or 


Pi + Pi 


P + Pr- 


(44) 


Pr <Pt ■ 
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Therefore 


pi+Pi= p + mm{p r ,p t }. (45) 

That is, if the the total stress-energy satisfies the NEC, then so does the perfect 
fluid portion. (Likewise, if the the total stress-energy violates the NEC then so does 
the perfect fluid portion.) 

Similarly, consider the DEC (dominant energy condition). The DEC requires 
(among other things) p ± p > 0. Now 

pi ~Pi = P~ ( Pr + Pt ^ - 1 \p r -p t \ = p - ma,x{p ri p t }; (46) 

Combining this with the previous result for the NEC, we see that if the DEC is 
satisfied for the total matter, then the DEC is satisfied for the fluid component. 

In contrast, consider the TEC (trace energy condition, p — 3p > 0). Note the 
TEC is not fundamental physics, and was for many years abandoned.^ But the 
TEC has undergone a recent resurgence of interest, due mainly to the fact that it 
can nevertheless usefully be used to characterize normal laboratory matter.^ We 
note 

3 {p r +pt) 1. , , . 

Pi — 3pt = P --- - \Pr-PtV (47) 

Then 

Pi — 3pf = p — (Pr + 2 Pt) +Pt~ ma x{p r ,p t } < p - (p r + 2 p t ). (48) 

So if the fluid component satisfies the TEC then the total matter satisfies the TEC. 
(In this particular case, the converse need not hold.) 

In summary, we see that the various classical energy conditions can quite easily 
be investigated in terms of the three-component model that we have set up to 
mimic any anisotropic fluid sphere. If there are any violations of the classical energy 
conditions, they are confined to the perfect fluid component of the model. 


6. Discussion 

So what have we accomplished here? On the one hand, we have constructed an 
explicit three-component model that is capable of mimicking any anisotropic fluid 
sphere, with two of the components (electro-magnetic field and scalar field) auto¬ 
matically being theoretically attractive and physically well behaved (eg, satisfying 
the classical energy conditions). It is only the perfect fluid component that has any 
risk of being “exotic” (ie, violating one or more of the classical energy conditions). 

On the one hand this is a very powerful result; on the other hand it is perhaps a 
little too powerful. Note that we are not claiming that any anisotropic fluid sphere 
is physically one of these three-component models, but are instead making the more 
modest claim that both the spacetime geometry and the total stress-energy can be 
successfully mimicked by one of these three-component models. 
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We feel that this is already an interesting observation, and that these three- 
component models may be of some theoretical (and maybe even practical) interest 
— certainly, since our construction is entirely generic, one can eliminate the need 
for making any ad hoc choices for the metric components or the equation of state, 
as is typically done for currently extant models of anisotropic fluid spheres. 
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Appendix A. General solution to the inversion problem 

The completely general solution to the inversion problem, now parameterized by an 
arbitrary function ft(r), is 


Pi = 
.n 2 _ 


Pr + Pt . 

2 ’ 

(A.l) 

ma x{p r — p t , 0} + ft 2 ; 

(A.2) 

ma x{p t — p r , 0} + ft. 2 ; 

(A.3) 

P~\ \Pr~Pt\ ~ h 2 . 

(A.4) 


Note that at the surface of the anisotropic fluid sphere we now have \E(r~)\ > 
|ft(r“)|. So if |ft(rj)| > 0, and if we are to embed this model into an exterior 
Schwarzschild geometry, then there must be a compensating shell of electric charge 
density at the surface. This is merely a feature of the model, it is in no sense 
problematic. 

With regards to the NEC, note that in particular 

Pr + Pt 1 | i,2 m r\ 

pi+Pi=p-\ -^- 2 \Pr~P t \~ h ■ ( A - 5 ) 

So in this case we have either 


Pi +Pi = P + Pt - h 2 ; 


Pr > Pt ■ 


(A.6) 
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or 

Pr < Pt • Pi + Pi = P+Pr - h 2 . (A.7) 

Therefore 

Pi +Pf = p + min{p r ,p t } - h 2 . (A.8) 

Note that statements about the NEC are now trickier: 

• Violation of the total NEC now implies violation of the NEC for the fluid 
component. 

• Violation of the fluid component NEC no longer implies violation of the 
total NEC. 

• Even if the total stress-energy satisfies the NEC, one can now always (by 
choosing h 2 large enough) force the fluid component to violate the NEC. 
(Or force W{ = —1, on the border of NEC violation. Note that a fluid 
component satisfying pf + Pi = 0 is now not cosmological constant because 
it is interacting with other forms of matter.) 

Regarding the DEC (the dominant energy condition), which requires pip > 0, 
let us consider 

Pi ~ Pi = P ~ ^ Pr + Pt ^ - ^\p r - Pt\ ~ h 2 = p - ma x{p r ,p t } - h 2 . (A.9) 

So if the DEC is satisfied for the total matter, one can always choose h to make the 
fluid “stiff”: pt = pf. 

When it comes to the TEC (the trace energy condition, p—3p > 0, see references 
l42l and S3), we have 

Pi — 3p f = p - - - - \p r -Pt\-h . (A.10) 

Then 

Pi - 3p { = p - (p r + 2 p t ) +p t - max{p r , p t } - h 2 < p - (p r + 2 p t )\ (A.11) 

So if the perfect fluid component satisfies the TEC, then the total matter content 
satisfies the TEC. 

Now consider the Strong Energy Condition (SEC), which (in addition to the 
NEC) states that p + 3p > 0. Then 

. q i 3(p r + Pt) 1 | | ,2 

Pi + opt = p~\ --- - \Pr ~ Pt\ ~ h 

. , . ■, . Pr+Pt 1, , ,2 

= P+(Pr +Pt) 4-2-2 \ Pr ~ Pt ' - h 

= P + ( Pr + Pt) + min {p r ,p t } - h 2 

= p + ( Pr + 2 p t ) -pt + min{p r ,p t } - h 2 

< p + (p r + 2pt). (A.12) 

So if the fluid component satisfies the SEC, then the total matter content satisfies 
the SEC. 
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